home *** CD-ROM | disk | FTP | other *** search
/ SGI Freeware 1999 August / SGI Freeware 1999 August.iso / dist / fw_plotutils.idb / usr / freeware / share / ode / population.ode.z / population.ode
Encoding:
Text File  |  1998-10-28  |  608 b   |  32 lines
  1. # This is an example from population biology (the predator-prey
  2. # equations of Lotka and Volterra).  You may run it by doing:
  3. #
  4. #    ode < population.ode | graph -T X -C
  5. # 
  6. # or alternatively, to get a real-time plot,
  7. #
  8. #    ode < population.ode | graph -T X -C -x 0 10 -y 0 3
  9. #
  10. # The plot shows the population of the prey as a function of time.
  11. # The curve oscillates, because as the prey die out, the predators
  12. # starve.
  13.  
  14. # The differential equations are:
  15. # x' = (A - By) x
  16. # y' = (Cx - D) y
  17. # A,B,C,D > 0
  18.  
  19. x' = (A - B*y) * x
  20. y' = (C*x - D) * y
  21.  
  22. A = 1
  23. B = 1
  24. C = 1
  25. D = 1
  26.  
  27. x = 3
  28. y = 1
  29.  
  30. print t, x
  31. step 0, 10
  32.